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We use cosmological constraints from current data sets and a figure of merit approach in order 
to probe any deviations from general relativity at cosmological scales. The figure of merit approach 
is used to study and compare the constraining power of various combinations of data sets on the 
modified gravity (MG) parameters. We use the recently refined HST-COSMOS weak-lensing tomog- 
raphy data, the ISW-galaxy cross correlations from 2MASS and SDSS luminous red galaxy surveys, 
the matter power spectrum from SDSS-DR7 (MPK), the WMAP7 temperature and polarization 
spectra, the baryon acoustic oscillations from Two-Degree Field and SDSS-DR7, and the Union2 
compilation of type la supernovae, in addition to other bounds from Hubble parameter measure- 
ments and big bang nucleosynthesis. We use three parametrizations of MG parameters that enter 
the perturbed field equations. In order to allow for variations of the parameters with the redshift 
and scale, the first two parametrizations use recently suggested functional forms while the third is 
based on binning methods. Using the first parametrization, we find that the CMB + ISW + WL 
combination provides the strongest constraints on the MG parameters followed by CMB + WL or 
CMB + MPK + ISW. Using the second parametrization or the binning methods, we find that the 
combination CMB + MPK + ISW consistently provides some of the strongest constraints. This 
shows that the constraints are parametrization dependent. We find that adding up current data sets 
does not improve consistently the uncertainties on MG parameters due to tensions between the best- 
fit MG parameters preferred by different data sets. Furthermore, some functional forms imposed 
by the parametrizations can lead to an exacerbation of these tensions. Next, unlike some studies 
that used the CFHTLS lensing data, we do not find any deviation from general relativity using 
the refined HST-COSMOS data, confirming previous claims in those studies that their result may 
have been due to some systematic effect. Finally, for all the parametrizations and binning methods 
used, we find that the values corresponding to general relativity are within the 95% confidence level 
contours for all data set combinations. 

PACS numbers: 95.36.+x,98.80.Es,98.62.Sb 

I. INTRODUCTION 

One of the routes to understanding the origin of cosmic acceleration is by testing whether the acceleration is due 
to a dark energy component pervading the universe or some modification to general relativity (GR) at cosmological 
scales of distances. The study of the growth rate of large scale structure can help to make such a distinction and 
has been the subject of a number of analyses in the literature. Some first papers (see for example [H-[l]) studied 
inconsistency between constraints on parameters from the growth rate and constraints from the expansion history 
to detect deviations from GR using observations. These were followed by another approach (see for example [El-l46j) 
that based their analysis on various growth parameters that take distinctive values for distinct gravity theories. For 
example, in GR, some parameters take the unity value, and deviation from it would indicate deviations from GR. It 
is worth mentioning that some useful pioneering work on the growth rate parametrization in for example (47l - [50j had 
been done for a purpose other than seeking to distinguish between dark energy and modified gravity (MG). Finally, it 
is informative to refer to all the early work about post-GR parametrizations that had been proposed to test deviations 
from GR at the level of solar system tests, for example, see [HJ and references therein. 

Recently, some studies used current available data sets in order to constrain the MG parameters and to see if there 
is an indication of any deviation from GR, but this has been in vain so far [53- 5fjj . 

Iln this paper, we calculate and compare constraints and figures of merit (FoM) on MG parameters based on 
func- tional forms or binning methods, using combined data sets of baryon acoustic oscillations (BAO) measurements 
from Two-Degree Field and DR7 Sloan Digital Sky (SDSS) surveys pJEU, the WMAP 7-year (WMAP7) cosmic 
microwave background (CMB) temperature (TT) and polarization (TE) spectra [62[, and the supernovae Union2 



* Electronic address: jnd041000@utdallas.edu 
' Electronic address: jam042100@utdallas.edu 
t Electronic address: mishak@utdallas.edu 



2 



compilation, which includes the 557 Type la SNe of the Supernovae Cosmology Project (SCP) (63|, and references 
of other compiled supernovae therein, the matter power spectrum (MPK) from SDSS DR7 [6(3], the integrated Sachs 
Wolfe (ISW)-galaxy cross correlations with the 2MASS and SDSS luminous red galaxy (LRG) surveys |64l. l65| . and 
the recently refined Hubble Space Telescope (HST) Cosmic Evolution Survey (COSMOS) weak- lensing tomography 
analysis in [66]. 



II. METHODOLOGY 



A. Unmodified Growth Equations 

For convenience, we will work in the conformal Newtonian gauge, where the perturbed Friedmann- Lemaitrc- 
Robertson- Walker metric is written as: 

ds 2 = a(r) 2 [-(l + 2^)dr 2 + (1 - 2^)dx i dx i }, (1) 

where a(r) is the scale factor normalized to one today, the x^s are the comoving coordinates, and r is conformal time. 
(f> and tp are scalar potentials describing the scalar mode of the metric perturbations. 

We can get two very useful equations that describe the evolution of the scalar potentials by using the first-order 
perturbed Einstein equations working in Fourier k-space. Combining the time-time and time-space equations gives 
the Poisson equation. Taking the traceless, space-space component of these equations lets us relate the two potentials. 
Respectively, we have: 

k 2 (j> = -4irGa 2 J2 P iAi (2) 

i 

k 2 {4>-<i>) = -12^Ga 2 ^ Pl (l + Wl )y, (3) 

i 

where pi and 7Tj are the density and the anisotropic stress, respectively, for matter species, i. Aj is the gauge- 
invariant, rest-frame overdensity for matter species, i, the evolution of which describes the growth of inhomogeneities. 
It is defined by: 

Ai = Si + 3n^, (4) 

where % = a /a is the Hubble factor in conformal time, and for species i, Si = 5pi/p is the fractional overdensity and % 
is the heat flux and is related to the divergence of the peculiar velocity, 9i , by 9i — 1+ ^ . Enforcing the conservation 
of energy-momentum on the perturbed matter fluid, these quantities for uncoupled fluid species or the mass-averaged 
quantities for all the fluids evolve as described in [67] : 

6 = -kq + 3{l + w)<j) + 3H(w~c 2 s )5 (5) 

| = -H(l 3w)l - T ±-l + c 2 s S-(l + w )l + { l + W )^. (6) 

Above, w — pj p is the equation of state and c 2 — Sp/Sp is the sound speed. Combining these two equations, we can 
express the evolution of A by: 



A =3(l + w)<f> + mwA - 



m-^— + k 2 + 3{n 2 -n 



Equations ^ , Q , © , and © are coupled to one another, combining them, along with the evolution equations for 
a(r), we can describe the growth history of structures in the universe. 



B. Modifications to the growth equations for detecting deviations from general relativity. 

Recently, a lot of attention has gone into detecting deviations from general relativity by parametrizing both mod- 
ifications to Poisson's equation, ([2]), as well as the ratio between the two metric potentials and c in the perturbed 
Friedmann-Lemaitre- Robertson- Walker metric (called gravitational slip by [Hj]), see f° r example [ill [531457T ]. In this 
paper, we use two continuous parametrizations and a binned parametrization that allow for modifications of Eqs. @ 
and ([3]) both directly and indirectly. 
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1. Examples of parametrization with functional forms for time and scale dependencies 

The first parametrization with a functional form that we use was proposed by 54]. It allows modifications to Eqs. 
(|2"|) and ^ that evolve monotonically in both time and scale, and makes no assumptions as to the time when a 
deviation from general relativity is allowed. These modifications are as follows: 



k 2 cj) = —AkGcl 2, piAi Q (8) 



e^-Rcj)) = -12nGa 2 J2pi(l + Wi )^-Q, (9) 

i 

where Q and R are the MG parameters. The parameter Q represents a modification to the Poisson equation, while 
the parameter R quantifies the gravitational slip (at late times, when anisotropic stress is negligible, R — ip/ft)- These 
parameters are parametrized to evolve in time and scale as: 

X(k,a) = \x Q e-^ k " + X tx (l ~e- k l k ") - lj a s + 1, (10) 

where X denotes either Q or R. So the model parameters which can be used to detect deviations from GR are now: 
Qo, Ro, Qoo, Roo, k c , and s. The parameters s and k c parametrize time and scale dependence respectively, with GR 
values s = and k c — oo. Qo and Rq are the present-day superhorizon values while Qoo and Roo are the present-day 
subhorizon values of the Q(k,a) and R(k,a) , all taking GR values of 1. 

The second parametrization we use comes from [57j and is used in their code MGCAMB [68j]. Here, a modification to 
Eq. @ is done indirectly by defining a modified field equation containing the parameter p, as well as, a gravitational 
slip parameter, r/. Explicitly these modifications are: 

fc 2 ?A = -ATiGa 2 ^p l A l p(k : a). (11) 

i 

7 = V(k,a) (12) 

V 

The parameters, {p, 77} are allowed to have a redshift dependence where they are fit to constant values below some 
transition redshift, z s , (note that this is not to be confused with the source galaxy distribution peak redshift) and 
then make a transition to a GR value of 1 as following a hyperbolic tangent function with a transition width, Az: 

= ^ f 1 + tanh i-ii) + m . (14) 



2 V Az 
Also, the parametrization uses one other parameter, S(/c, a), defined by: 

but with only a redshift dependence here in this subsection. As we discuss in Sec. II C below, this parameter is useful 
to break degeneracies between the other parameters. 



2. Example of parametrization with binning in time and scale 

Third, in order to allow for simultaneous scale (wavenumber) and redshift (time) dependences in the {^{k, a), a)) 
parametrization above, we use variants of a 2 x 2 binned approach (called pixillation in 57]). We use two bins for 
scale dependence as 0.0 < k < k x and k x < k < 00, and where we use an analysis with k x = 0.01 hMpc -1 and the 
other with k x = 0.1 hMpc -1 . In both cases, we bin for redshift dependence as < z < 1 and 1 < z < 2. We also use 
a third redshift binning method with < z < 1.5 and 1.5 < z < 3 and the scale bin separator at k x = 0.01 hMpc -1 
for the data sets in this paper. We use combinations of the parameters pi and with i = 1...4 for each [k, z] bins in 
a 2 x 2 grid as shown in Table fl] Transitions between the redshift bins behave as a hyperbolic tangent function with 
a transition width of Az = 0.05 (almost a step), while transitions between scale bins are straightforward steps. 
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Redshift bins 


Scale bins 


0.0 < z < 1, 1.5 


1, 1.5 < z < 2, 3 


0.0 < k < k x 


Si 


^2, S 2 


k x < k < oo 


M3, S 3 


/i4, E4 



TABLE I: The layout of the binned parametrizations. Specifically, for the first two binned methods this involves using {pi, £1} 
for the < z < 1 and 0.0 < k < k x bin, {^2, £2} for the 1 < 2 < 2 and 0.0 < k < k x bin, {^3, S3} for the < z < 1 and 
k x < k < 00 bin, and {/m, E4} for the 1 < z < 2 and k x < k < 00 bin, and the third binned method uses {/j,i, £1} for the 
< z < 1.5 and 0.0 < fc < fa bin, {^2, £2} for the 1.5 < z < 3 and 0.0 < fc < fa bin, {^3, £3} for the < z < 1.5 and 
fa < fc < 00 bin, and {^4, £4} for the 1.5 < z < 3 and k x < k < 00 bin. 



C. Relations between parametrizations and degeneracies among parameters 

We now turn our discussion to relating the different parametrizations to one another. We will focus on the {Q, R} 
and {/x, r), £} parametrizations above, as well as, a third parametrization using {Q, V} used in|56|. When relating the 
MGCAMB and {Q, V} parametrizations to that of the first, we give the relations during matter domination, assuming 
zero anisotropic stress. 

. = | = 2^-1 (16) 

m = Q = 2Q-V. (17) 

We use these relations later in order to infer constraints on some parameters from the ones obtained from the 
constraints derived from the data sets. 

As previously pointed out (54l - [57| some of these parameters are degenerate and it is worth introducing parameters 
that alleviate these degeneracies. First, the parameters Qo,oo and Rq.oo are degenerate along a direction that we call 
-Do, 00 = Qo,oo(l + Rq.oo) /2 this degeneracy is shown with a curve in Fig. [TJ We proceed by varying Qo,oo and Do, 00 
and infer values of Rq i00 from the two parameters. Similarly, fi and 77 are degenerate along curves of £ — constant so 
we vary /i and S to get rid of this degeneracy. Using the parameters D and £ is also useful because observations of 
the weak-lensing and ISW are sensitive to the sum of the metric potentials 4> + tp and its time derivative respectively. 
Thus observations are able to give us direct measurements of these parameters. 



H = QR = V 



D. Figures of merit 



We calculate the figures of merit for the parameter pairs used for each parametrization considered. The FoM 
approach has been used previously, for example, to quantify the constrai ning po wer of cosmological data on the 
dark-energy equation-of-state parameters, using simulated future data [ill [69l - l73j . or current data [74T - [76l |. For a 
pair of parameters, the FoM can be defined to be proportional to the inverse of the 95% confidence limit area. For 
more parameters, this can be generalized to be proportional to the inverse of the volumes and supervolumes of the 
(super)ellipsoids, again representing the confidence limits in these parameter hyperspaces. A variety of choices for 
the constants of proportionality has been made in previous work and we simply chose to set it to one. Assuming a 
Gaussian distribution, the FoM for two parameters is inversely proportional to the area of the ellipse comprising the 
contour of constraints for the two parameters and so we define our FoM [74j , 

FoM = {detC)- 1/2 (18) 

where C is the 2D covariance matrix for the two parameters • This is similar to the FoM used by the Dark Energy 
Task Force for the equation of state parameters. Obviously, tighter constraints lead to higher FoMs. 

For two uncorrelated parameters with Gaussian uncertainties, the FoM does not provide any more information than 
the individual uncertainties in describing the constraining power of the data sets [69j, but our results are not purely 
Gaussian, and our parameters for which we show FoM's are not always uncorrelated. The FoM's are very informative 
to tell us which data set is most significant in constraining the parameters and in our study here, were useful in 
indicating some tension between the data sets. It is of course understood that the FoM as we have defined it here, 
is only an approximation with non-Gaussian uncorrelated parameters, but accurate within ~ 10% in the worst cases 
[z3|. The FoM is less accurate for two degenerate parameters, but still provides a good approximation to compare 
the constraining power of data sets. In fact, at first glance, when looking at the constraints in the tables and figures 
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0.5- 
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CMB+MPK+WL 

CMB+MPK+ISW 

- - CMB+MPK+ISW+WL 




-CMB+MPK+WL 
-CMB+MPK+ISW 
CMB+MPK+ISW+WL 




0.5 
V..u, 



CMB+MPK+ISW+WL, k = 0.01 
CMB+MPK+ISW+WL, k . «. 




- CMB+MPK+ISW+WL, k c = 0.01 
CMB+MPK+ISW+WL, k = ~ 




FIG. 1: TOP: 68% and 95% C.L. for the constraints on Qo and R . TOP-LEFT: Constraints with Q and R scale independent 
(k c — oo). To show how the parameter Do actually represents the degeneracy direction of these parameters, the solid black line 
is the line Do = 1. TOP-RIGHT: A comparison of constraints when Q and R are allowed either scale dependence (k c = 0.01) 
or scale independence using CMB, ISW, MPK, and WL. BOTTOM-LEFT: Inferred 68% and 95% C.L. constraints for the 
parameters Q or Ei and V or fii in the second and third parametrizations 56, 57]. These constraints are inferred by using 
constraints on the parameters Qo, Ro, and Do from the first parametrization [5f|. For all of these constraints we set k c — oo 
implying scale independence of the parametrization (|10|) . BOTTOM- RIGHT: We compare inferred constraints on Q or Ei and 
V or fi\ using CMB, ISW, MPK, and WL when scale dependence is allowed in the first parametrization to constraints when 
the parametrization is scale independent. 



below, without the FoM, it is more difficult to determine which combination of data sets is more sensitive because at 
times there is a shift in the allowed parameter space for the different combinations of data sets. 

This metric is useful in comparing the constraining power of various combinations of data sets that we used, 
namely WMAP7 temperature and polarization spectrum (CMB), the matter power spectrum (MPK), ISW-galaxy 
cross correlations (ISW) and HST COSMOS weak-lensing tomography (WL). In each case, we use Union2 supernovae 
(SN), BAO, BBN, AGE and H , as well. 



III. DATA SETS FOR OBSERVATIONAL CONSTRAINTS 

We use Type la SNe from the Union2 compilation (6^. BAO from Two-Degree Field and SDSS DR7 H [pL 
and WMAP7 spectrum. The WMAP7 data, MPK from SDSS DR7 [H, the ISW-galaxy cross correlations [gl i|| 
and COSMOS weak-lensing tomography [66[ will constrain the growth as described below. We compare to the 
full temperature (TT) and polarization(TE) spectra using the WMAP 7-year data release using the likelihood 
routine provided by the WMAP team (62|. Next, for the BAO, we follow the most recent comparison IMl and 
define the ratio of the sound horizon, rgjzd) to the effective distance, Dy as a fit for SDSS. We also add the prior 
H a = 74.2 ± 3.6 km/s/Mpc given by [78], and the priors on the age of the Universe (AGE) 10 Gyrs<AGE< 20Gyrs, 
and fifth 2 — 0.022 ± 0.002 from big bang nucleosynthesis (BBN) . A possible deviation from GR is allowed by the 
modified Poisson equation and inequality between <fi and ip. The growth of structure is important in constraining the 
model parameters for a deviation from GR in this modified description. For the total matter power spectrum today, 
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FIG. 2: We compare the 68% and 95% C.L. constraints on the Qo and Ro (LEFT) as well as the inferred parameters Q or Ei 
and V or m (RIGHT) for CMB+MPK (solid line) and CMB+WL (dashed line) and CMB+MPK+WL (dash-dotted) to show 
the tension between these data sets. Notice the combined constraints are outside those of the individual data sets. 



for example, see [EiJ, we use 

P(fc) = 27 r 2 (^) 3 AL, (19) 

where A tot = (p c A c + pf,A.f,)/{p c + pt) is the gauge-invariant, rest-frame, overdensity for cold dark matter and baryons 
combined. Recall, that the matter power spectrum is sensitive to the growth of structure by constraining p through 
A [54], [HI, WMAP7 data also provides large scale ISW measurements. The ISW effect describes the evolution 
in time of the gravitational potentials. The photons from the CMB gain energy falling into varying potential wells, 
and may not return it all on exit. This causes a secondary anisotropy of the CMB and can be analyzed directly with 
the derivatives of 4> an d ip to give an indication of deviation from GR, as described in [H3, It is difficult 

to distinguish these CMB temperature anisotropics signatures from rest of CMB anisotropics because it has 1/10 
the amplitude, so we need to correlate this effect to the large scale structure of the universe. For this we use the 
ISW-galaxy cross correlations which include the use of 2MASS and SDSS LRG galaxy autocorrelations explored in 
[H, [65[ and recently used by 0, HH . 

We further constrain the parameters with expansion history and the effect of the growth of structure from the weak- 
lensing tomography shear-shear cross correlations of the HST COSMOS survey recently compiled by Schrabback et. 
al. 16611 . The weak-lensing constrains the growth of structure through the combination of <j> and tp as described in 
[E3, 13 H3] ■ The authors of [IB] perform a refined analysis of the HST COSMOS survey of in combination with 
the COSMOS-30 photometric redshift catalogue provided by [soj . The shear-shear cross correlations were calculated 
between 6 redshift bins, 0.0 < z < 0.6, 0.6 < z < 1.0, 1.0 < z < 1.3, 1.3 < z < 2.0, 2.0 < z < 4.0 and a sixth 
bin that contains all faint galaxies with a numerically estimated redshift distribution from 0.0 < z < 5.0, see Fig. 6 
of [66| , Certain exclusions were made as described in (66|, such as the LRGs to avoid G-I intrinsic alignment bias 
and the lowest angular theta bin because of model uncertainties. Only bright galaxies (i < 24) were used in the first 
bin (z < 0.6). Also, autocorrelations were not used in bins 1 — 5 to reduce the effect of I-I intrinsic alignments. For 
comparing the shear-shear cross correlations, we modify the publicly available code for the COSMOS 3D weak-lensing 
built by Lesgourgues et. al. UK to incorporate the shear cross correlations as seen in, for example, (6(| (see endnote 

For the ISW-galaxy cross correlations we modified our code as described in previous works [H, [5(1 H3| and refer the 
readers to those works for an extended discussion of the implications of these modifications. Here, however, we use 
the newly refined HST-COSMOS 3D weak-lensing tomography by [66] and describe the modeling of this data below 
(We note that during the preparation of this manuscript, [59( released a paper fitting MG parameters using this data 
set as well). 

In order to model the weak-lensing, the shear cross correlation functions & _(#) between bins k,l are given by: 

1 f°° 

£+,_(*) = ^J dl i J 0A (£9)P* l (£), (20) 

where J„ is the n*' l -order Bessel function of the first kind, I is the modulus of the two-dimensional wave vector, and 
P^ 1 is the convergence cross-power spectra between bins k, I. These cross-power spectra are computed from integrating 
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the three-dimensional (nonlinear) power spectrum P$ 



vm—srl ^ p 4 W * (21) 

with comoving radial distance x, comoving distance to the horizon Xh, an d comoving angular diameter distance Jk{x)- 
We note here that the effect of the MG parameters on the input power spectrum from the transfer function has been 
implemented in our weak-lensing code, while the effect of the growth has been accounted for by our modifications to 
CAMB [83j], and both are included in P$ of Eq. (|21l) . Also, following (66|, we weight the geometric lens- efficiency factors 

g k (x) = r d X %( X ') fK ^'~ X) , (22) 

corresponding to the redshift distributions pk for the two compared redshift bins. We use the 160 x 160 dimensional 
covariance matrix provided by Schrabback et. al. [66j j and for the corresponding weak shear correlations. We also 
apply the correction to the inverse covariance, Cod -1 , seen in [gl, HU, as 

Cov - 1 = 0.4390 Gov" 1 . (23) 

for 288 independent realizations and a 160 dimensional data vector. Following instruction from (66|, we weight the 
predictions for the k th 9 bin, 9^, with logarithmic spaced upper and lower limits 9k^ max and 9k, m in, respectively, based 
on the number of galaxy pairs, ./V for a given 9 as 

rfl»,«« N(e')Z + -(9')d9' 
£(0 fc ) = Jflfc — (24) 

J*::::: 

with N(9) oc 0(0.0004664 + 9(0. 0044118 - 8.90878 x 10 -5 6>)) where the units of 9 is arcminutes. jMIH- Finally^ the 
10% uncertainty in the numerical estimate of the galaxy redshift distribution of bin 6 is applied, following [66j, as 
P&{ z i fz) = P6(fzz), where f z is a nuisance parameter marginalized in the results below. 



IV. RESULTS AND DISCUSSION 



Monte Carlo Markov chains (MCMC's) are used to compute the likelihoods for the parameters in the model. This 
method randomly chooses values for the above parameters and, based on the x 2 obtained, either accepts or rejects the 
set of parameters via the Metropolis- Hastings algorithm. When a set of parameters is accepted it is added to the chain 
and forms a new starting point for the next step. The process is repeated until the specified convergence is reached 
(via the Raftery-Lewis convergence criterion). We perform the best-fit of the parameters in the various models by 
X 2 minimization, i.e. minimizing for example x 2 = Xsn + X 2 B ao + Xcmb + xjsw + Xwl + Xmpk via a maximum 
likelihood analysis. To get best-fits for the MG parameters of a given parametrization we use, as discussed above 
modified versions of the publicly available ISW and Weak-lensing Likelihood code (64j, HST-COSMOS weak-lensing 
code [81j . CAMB (83|, and CosmoMC [84j . For these modified codes, see [s3 | . We fix the dark energy equation of state 
w = — 1 assuming a ACDM expansion. Then in addition to varying the MG parameters for a given parametrization, 
we vary the six core cosmological parameters: ft^h 2 and the f2 c /i 2 , the baryon and cold-dark matter physical density 
parameters, respectively; 9, the ratio of the sound horizon to the angular diameter distance of the surface of last 
scattering; T re i, the reionization optical depth; n s , the spectral index; and lnl0 lo A s , the amplitude of the primordial 
power spectrum. These parameters are allowed to vary within the wide range of values that are defaulted in the 
params.ini file of CosmoMC. Additionally, as discussed in the previous section, when using the WL data set we 
marginalize over the nuisance parameter f z . We provide various 95% confidence limits on the MG parameters along 
with the corresponding contour plots. We also calculate the figure of merit (FoM) for the parameter combinations for 
the parametrization considered. 

As described above, for the first parametrization, i.e. Eqs. © and ©, we follow the technique of varying Qo,oo and 
-Do.oo = Qo,oo (1 + Ro.oo) /2, while inferring values of i?o,oo from the other two parameters. We allow these parameters 
to vary < Qo,oo < 10, < -Dq.oo < 10, additionally ruling out parameter combinations where i?o,oo < — 1- We fit 
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Constraints for the parameters {Qo.oo, Ro,oo, Do.oo} 


k c — oo 


data set 


FOM Qo;flo 


FOM Qo;Do 


Qo 


Ro 
(derived) 


Do 




Roo 

(derived) 


D x 


i\ r 1 1 ft r ~r~i t /~ 

CMB, MPK 


5.53 


8.81 


[0.69, 2.79J 


[—0.34, 2.04] 


[0.79, 1.62] 








CMB, ISW 


6.15 


9.80 


[0.70, 2.76] 


[-0.33, 1.85] 


[0.79, 1.55] 








CMB, WL 


6.44 


12.15 


[0.57, 2.54] 


[-0.30, 2.24] 


[0.76, 1.47] 








CMB, ISW, WL 


8.46 


14.96 


[0.67, 2.53] 


[-0.33, 1.84] 


[0.75, 1.34] 








CMB, MPK, ISW 


7.53 


11.97 


[0.80, 2.73] 


[-0.35, 1.62] 


[0.76, 1.46] 








CMB, MPK, WL 


5.08 


7.96 


[0.73, 2.80] 


[-0.34, 2.13] 


[0.82, 1.65] 








CMB, MPK, ISW, WL 


7.09 


10.06 


[0.83, 2.82] 


[-0.37, 1.60] 


[0.77, 1.56] 








k c = 0.01 


data set 


FOM Qo;Ro 


FOM Qo;Do 


Qo 


-Ro 
(derived) 


Do 


QoG 


Roo 

(derived) 


Doo 


CMB, MPK, ISW, WL 


4.07 


4.53 


[0.57, 3.35] 


[-0.79, 1.67] 


[0.29, 1.60] 


[0.46, 4.01] 


[-0.82, 2.28] 


[0.26, 2.00] 





Inferred parameter constraints on the parameters Q, Ei and V, from parameters Qo, Ro, and Do 


k c — oo 


data set 


FOM 5; v 


0,Ei 


V,Mi 


CMB, MPK 


8.81 


[0.79, 1.62] 


[-0.90, 1.79] 


CMB, ISW 


9.80 


[0.79, 1.55] 


[-0.86, 1.57] 


CMB, WL 


12.15 


[0.76, 1.47] 


[-0.72, 1.60] 


CMB, ISW, WL 


14.95 


[0.75, 1.34] 


[-0.78, 1.38] 


CMB, MPK, ISW 


11.97 


[0.76, 1.46] 


[-0.89, 1.42] 


CMB, MPK, WL 


7.96 


[0.82, 1.65] 


[-0.90, 1.91] 


CMB, MPK, ISW, WL 


10.06 


[0.77, 1.56] 


[-0.98, 1.57] 


k c = 0.01 


data set 


FOM S;V 




V,Mi 


CMB, MPK, ISW, WL 


4.53 


[0.29, 1.60] 


[-2.46, 1.48] 



TABLE II: TOP: 95% C.L. for the parameters Qo, Ro,Do, Qoo, Roo, and Doo in the first parametrization [541 ]. For these 
constraints, we marginalize over the parameter s which characterizes the redshift (time) dependence of the parametrization. 
BOTTOM: Inferred 95% C.L. constraints for the parameters Q or Ei and V or fii in the second and third parametrizations 

5(1 H3- These constraints are inferred by using constraints on the parameters Qo, Ro, and Do from the first parametrization 

11. 



both scale- independent and scale-dependent forms of this parametrization, in both cases we allow for time dependence 
by marginalizing over the parameter s as described earlier. For scale-independent modifications, we fix k c = oo and 
of course vary only the parameters Qo a nd Do, while for scale-dependent modifications we fix k c = 0.01, allowing all 
four parameters Qo.oo and -Do.oo to vary. The 95% C.L.'s on these parameters and corresponding FoM's for various 
combinations of cosmological probes are given in Table [TTJ We also plot the 68% and 95% 2D contours for Qo and Ro 
in Fig. [T] We find in all cases consistency with general relativity at the 95% confidence level. 

We note here that we do not find the slight 2a deviation from GR on the parameter Qo seen in [54| when combining 
all the data sets. Rather in all cases for this parametrization we find parameter values consistent with general relativity 
at the 95% confidence level and a little looser than the results given in 54j . We were however able to reproduce almost 
exactly the same results of [54| but only when fixing the core cosmological parameters to their WMAP7 [62| best-fit 
values and allowing only the MG parameters to vary. We believe that the method of varying all the cosmological 
parameters simultaneously with the MG parameters gives more realistic and conservative results and those are the 
results that we reported in our Table ITT1 and Figs. Q]and[2] 

Next, using the constraints we attained on the parameters of the first parametrization [HiJ we can, with the help 
of Eqs. (fTo) and (fTT)) . infer constraints on the low redshift large scale parameters (^i and Ei below) and {Q and V 
in the low redshift and small k— bin therein). We list these inferred 95% CL's and their FoM's in the second part of 
Table ITU and plot the 68% and 95% 2D contours in Fig. [TJ We verified that the shapes of the contours plotted for the 
inferred parameters reflect the physical priors placed on the noninferred parameters as discussed above. Again, these 
results show consistency with general relativity at the 95% confidence level. 

Using the first parametrization, we find that the combination CMB+ISW+WL provides the strongest constraints 
on the parameters. We also observed that there seems some tension between current WL and MPK data. This 
reveals itself particularly well when comparing the FoM's of fits that include both these data sets to those that 
include only one of them eg. CMB+MPK+WL vs CMB+WL. Similarly, the combinations CMB+ISW+WL or 
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Constraints on /io 5 ??o? and Eo for redshift dependence 








Z s = 1 








Zs = 2 










Az 


= 0.05 










data set 


FoM 


Mo 




Vo 

(derived) 


FoM 


Mo 


So 


Vo 

(derived) 


CMB, WL 


94.61 


[0.621, 1.622] 


[0.948, 1.100] 


[0.233, 2.328] 


162.5 


[0.645, 1.351] 


[0.946, 1.080] 


[0.467, 2.206] 




i nn n 

100.9 


[0.610, l.otwj 


rnnKK i i no! 
[0.900, 1.103J 


[0.29/, 2.350J 


1 / i\0 


rn a i o i onQl 

[0.618, i.zyoj 


[0.900, 1.0/9J 


[0.O42, l.6l(\ 


CMB, ISW 


114.2 


[0.590, 1.427] 


[0.954, 1.105] 


[0.415, 2.505] 


196.0 


[0.635, 1.228] 


[0.952, 1.082] 


[0.637, 2.251] 


CMB, MPK, WL 


104.5 


[0.624, 1.549] 


[0.960, 1.108] 


[0.317, 2.334] 


179.4 


[0.643, 1.331] 


[0.959, 1.082] 


[0.514, 2.221] 


CMB, ISW, WL 


123.9 


[0.589, 1.399] 


[0.954, 1.098] 


[0.438, 2.499] 


209.2 


[0.627, 1.198] 


[0.954, 1.081] 


[0.658, 2.297] 


CMB, MPK, ISW 


132.9 


[0.594, 1.355] 


[0.956, 1.098] 


[0.487, 2.472] 


230.2 


[0.608, 1.157] 


[0.957, 1.076] 


[0.723, 2.372] 


CMB, MPK, ISW, WL 


132.5 


[0.592, 1.345] 


[0.961, 1.106] 


[0.520, 2.483] 


229.1 


[0.626, 1.182] 


[0.962, 1.081] 


[0.706, 2.307] 


Marginalizing over Az 


data set 


FoM 


Mo 


So 


Vo 

(derived) 


FoM 


Mo 


So 


Vo 

(derived) 


CMB, MPK, ISW, WL 


74.29 


[0.618, 1.371] 


[0.963, 1.226] 


[0.527, 2.471] 


194.0 


[0.674, 1.176] 


[0.957, 1.113] 


[0.706, 2.089] 




z s = 3 


z s = 0.76 








Az 


= 0.05 










data set 


FoM 


Mo 


So 


Vo 


FoM 


Mo 


So 


Vo 


CMB, WL 


223.9 


[0.709, 1.249] 


[0.944, 1.075] 


[0.578, 1.906] 


75.74 


[0.631, 1.676] 


[0.944, 1.119] 


[0.209, 2.276] 


CMB, MPK 


231.4 


[0.678, 1.214] 


[0.955, 1.076] 


[0.637, 2.063] 


83.27 


[0.629, 1.622] 


[0.952, 1.122] 


[0.254, 2.281] 


CMB, ISW 


250.9 


[0.724, 1.179] 


[0.949, 1.082] 


[0.684, 1.861] 


89.01 


[0.606, 1.531] 


[0.949, 1.121] 


[0.331, 2.410] 


CMB, MPK, WL 


225.8 


[0.709, 1.246] 


[0.952, 1.077] 


[0.602, 1.911] 


83.35 


[0.629, 1.625] 


[0.955, 1.124] 


[0.258, 2.302] 


CMB, ISW, WL 


288.1 


[0.717, 1.146] 


[0.951, 1.076] 


[0.721, 1.881] 


87.78 


[0.599, 1.539] 


[0.944, 1.115] 


[0.313, 2.424] 


CMB, MPK, ISW 


296.6 


[0.681, 1.109] 


[0.954, 1.077] 


[0.797, 2.045] 


96.48 


[0.600, 1.488] 


[0.948, 1.115] 


[0.358, 2.428] 


CMB, MPK, ISW, WL 


316.9 


[0.703, 1.127] 


[0.961, 1.076] 


[0.769, 1.961] 


95.59 


[0.603, 1.492] 


[0.955, 1.121] 


[0.368, 2.423] 


Marginalizing over Az 


data set 


FoM 


Mo 


So 


Vo 


FoM 


Mo 


So 


Vo 


CMB, MPK, ISW, WL 


272.8 


[0.729, 1.139] 


[0.956, 1.095] 


[0.758, 1.859] 


43.15 


[0.589, 1.559] 


[0.948, 1.298] 


[0.361, 2.689] 



TABLE III: The FoM for each data set and {/io, So} is given. 95% C.L. for the parameters /io, Vo, and So from the second 
parametrization [57, 68] for redshift dependence z 3 = 1, 2, 3 for different combinations of data sets. We also add the physically 
motivated transition redshift (from deceleration to acceleration in an Lambda cold dark matter model) given by 1 + zu-ans = 
(2n A /n,„) 1/3 « 1.76 (for example see Eq. (14) in [13) for comparison. The data used is described in the text above and 
indicated as: WMAP7 temperature and polarization spectrum (CMB), the matter power spectrum (MPK), ISW-galaxy cross 
correlations (ISW), weak-lensing tomography (WL). In all cases, the data used is combined with SN, BAO, BBN, AGE and 
H . 

CMB+MPK+ISW do better than the combination CMB+MPK+ISW+WL. In order to investigate this, we compared 
the best-fit models preferred by each combination of data sets finding that significant tensions exist between the best- 
fit values for the MG parameters when using this parametrization. There are no significant tensions between the 
core cosmological parameters. This tension is illustrated in Fig. [5] for both constraints on Qo and Rq as well as the 
inferred parameters Q or Ei and V or Mi- Also, in searching for the origin of this tension, we have both fixed the six 
core cosmological parameters to the WMAP7 values and additionally the time dependence parameter to s = 3, yet 
it persists, significantly. We also find that the second parametrization below does not exhibit this particular tension 
between WL and MPK but shows other tension although substantially less pronounced. 

For the second parametrization, i.e. the one used in MGCAMB 57], we derive uncertainties and FoM for n(k, a) 
and £(fc, a) to avoid the strong degeneracy between fi(k, a) and rj(k, a). In order to allow for the parameters to be 
redshift (time) dependent and scale (wavenumber) dependent, we use a binned approach to look for deviation from 
GR, with no functional form assumed for the parameters. First, however, we allow only redshift dependence in order 
to compare to the first parametrization [5_4J and look for deviation at late times, i.e. z < z s with a threshold redshift, 
z s = 0.76, 1, 2, or 3 and assume that at early times fj,(k, a) = 7y(/c, a) — I. As a remainder, the evolution in time 
is described by (TT5]) and (fl4"|) . We use the transition width in (IT31 and (fl"4"|) to be both fixed, Az = 0.05, so that 
the transition is almost a step function, as well as, allowed to vary in the interval [0.01, 0.5], becoming a nuisance 
parameter marginalized over. We allow these parameters to vary as follows: — 10 < Mo < 10, < So < 10, and ruling 
out parameter combinations where vo < — 1- For each data set used, our results are given in Table ILnl showing the 
FoMs and 95% confidence limits on /io and So, as well as the inferred constraints on r]o, for this case (i.e. redshift 
(time) dependence only). See Fig. [3] for comparisons of the 95% confidence limit contours. We find that the general 
relativity values, (/io — 1, Sq = 1), are within 95% confidence level for all data sets. The constraints from z s = 2 are 
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FIG. 3: TOP-LEFT:68% and 95% C.L. for the constraints on fio and So of the second parametrization [57], [68j with threshold 
redshift, z s = 1 and transition width fixed, Az = 0.05. All contours are fit for CMB, SN, BAO, BBN, AGE and H - Contours 
enclosed by a solid line include WL, contours enclosed by a dashed line include WL and MPK, and contours enclosed by a 
dash-dotted line include WL, MPK, and ISW. TOP-CENTER: Same as LEFT, but for threshold redshift z s = 2. TOP-RIGHT: 
Same as LEFT, but for threshold redshift z s = 3. BOTTOM-LEFT:68% and 95% C.L. for the constraints on fi and E of the 
second parametrization WH . [HI with threshold redshift, z s = 1 and transition width fixed, Az — 0.05. All contours are fit for 
CMB, MPK, SN, BAO, BBN, AGE and Ho. Contours enclosed by a solid line include WL, contours enclosed by a dashed line 
include ISW, and contours enclosed by a dash-dotted line include WL and ISW. BOTTOM-CENTER: Same as LEFT, but for 
threshold redshift z a = 2. BOTTOM-RIGHT: Same as LEFT, but for threshold redshift z s = 3. The point (1, 1) indicates the 
GR values. 



tighter than z s = 1 because the earlier threshold redshift allows more change in gravitational potential and growth 
from the data sets, also noted in [53] ■ We also add constraints where we set the threshold redshift z s — 3 taking 
advantage of the higher redshift (i.e. 2 < z < 3) WL data in HST-COSMOS. We know that the relevant redshift for 
cosmic acceleration is much smaller, but it is of interest to simply study the constraints obtained on the growth of 
structure using various available ranges of data. From the physically motivated transition redshift (from deceleration 
to accelerati on i n an Lambda cold dark matter model) given by 1 + Ztrans 

= {2n A /n m )^ 3 w 1.76 (for example, see 



Eq. (14) in [87jj with Q\ = 0.73, Cl m — 0.27), we add z s — 0.76 for comparison, also. The physically motivated 



transition redshift finds the weakest constraints and lowest FoM's due to less data below z s = 0.76, but all results are 
consistent with GR. 

The strongest constraints and highest FoM's we find come from using the combination CMB+MPK+ISW for 
both z s — 1, 2, whereas, for z s = 3 the tightest constraint comes from using all the data sets including lensing, 
i.e. CMB+MPK+ISW+WL because with z s = 3, we allow more WL data below the threshold redshift to constrain 
the MG parameters which also yields a longer constant value at late times. We find that the FoM for CMB+ISW 
is better than the one for CMB+MPK, which in turn is better than the FoM for CMB+WL. This shows that the 
current ISW data provides indeed the strongest constraints on the MG parameters. Also, as shown in Table [TTTT the 
allowance of the transition width, Az, to vary (i.e. [0.01, 0.5]) loosens the constraints and reduced the FoM for all z s 
models. Now, as found for the first parametrization, here also, there is a slight decrease in the FoM and weakening 
of the constraints (for z s = 1, 2) when WL is added to CMB+MPK+ISW, again, due to some tensions between the 
preferred MG parameters by different data sets. To investigate these tensions, here we also fixed the core cosmological 
parameters and left the MG parameters to vary. As a result, the FoM's increase consistently (but only moderately 
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FIG. 4: 68% and 95% C.L. for the parameters fj,i and Ei from the second parametrization [H?], EH) for redshift and scale 
dependence as a binned 2x2 parametrization with < z < 1, 1 < 2 < 2 and scale separation at k x = 0.01 h Mpc~ x . All 
contours are fit for CMB, MPK, SN, BAO, BBN, AGE and Ho- Contours enclosed by a solid line include WL, contours enclosed 
by a dashed line include ISW, and contours enclosed by a dash-dotted line include WL and ISW. The point (1, 1) indicates the 
GR values. 



in some cases) when adding a data set, indicating that some of the tension was reduced but not removed. In fact, a 
closer look at the best-fit MG parameters preferred by different data sets shows that, while reduced, a non-negligible 
scatter is still present. This explains the only moderate improvements obtained in some cases. This indicates that 
there are non-negligible tensions intrinsic to the MG parameters preferred by different data sets. 

Next, we calculate constraints and FoM's while we allow for simultaneous scale (wavenumber) and redshift (time) 
dependence by using variants of a 2 x 2 binned approach (called pixilation in see earlier discussion and Table 

Ufor details. We alternate between ISW and WL and their combination, while using CMB, MPK, BAO, SN, BBN, 
AGE and Hq in all cases. We provide the FoM's and 95% confidence limits on \ii and E^, as well as, the inferred ones 
on r\i in Table HVl for scale-dependence separators at k x = 0.01 h Mpc^ 1 and k x — 0.1 hMpc -1 , see Figs. H]and[5]for 
the comparisons in their C.L. contours. General relativity values for pi — 1 and E,; = 1 are within all constraints for 
Hi and Ei for all the combination of data sets. Unlike the work that uses a binned parametrization and CFHTLS 
[Hoi we do not find a deviation from GR with E 3 . But as expressed there as well, this is due to the fact that 
we are using the refined HST-COSMOS data by 66j which does not have the known systematic effect (i.e. residual 
systematic 'bump') as discussed in the papers using CFHTLS, 56, 57, 88, 89|. 



Another observation that can be made is that with the data sets used in this paper, the binning with scale separator 
at k x = 0.01 hMpc~ x better constrains {//i, E^} by the balancing of the number of data points within the bins, so 
within each bin, we are taking maximum advantage of the available data. Whereas, the wave number separator at 
k x = 0.1 h Mpc" 1 does not equally balance the available data, i.e. {p3, E3} and {/Z4, E4} are not well constrained 
by ISW because most ISW data points are below k x = 0.1 hMpc^ 1 , which allows WL to dominate the {ps, E3} and 
{/14, E4} bins. However, for bins with a scale separator at k x = 0.01 hMpc -1 , the FoM increases for {/i^, Ej} in each 
bin beginning with the lowest for WL, then ISW, and the highest for the combination of ISW and WL because of the 
balance of data points, see Table [TV] see Fig. H (again CMB, MPK, BAO, SN, BBN, AGE and H Q are included each 
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Constraints for {fii, E^} binned parametrization < z < 1, 1 < z < 2 





k x = 0.01 


k x = 0.1 


data sets 


WL 


ISW 


ISW, WL 


WL 


ISW 


ISW, WL 


FoMi 


75.92 


101.7 


101.2 


80.63 


99.47 


97.62 


Mi 


[0.559, 1.567] 


[0.480, 1.227] 


[0.494, 1.309] 


[0.587, 1.540] 


[0.542, 1.356] 


[0.535, 1.357] 


Ei 


[0.940, 1.138] 


[0.957, 1.145] 


[0.955, 1.139] 


[0.946, 1.141] 


[0.949, 1.136] 


[0.951, 1.139] 


(derived) 


[0.281, 2.815] 


[0.605, 3.486] 


[0.563, 3.347] 


[0.322, 2.644] 


[0.493, 2.939] 


[0.496, 3.004] 


FoM 2 


148.9 


176.3 


184.2 


161.7 


172.4 


180.5 




[0.614, 1.456] 


[0.563, 1.289] 


[0.580, 1.294] 


[0.633, 1.413] 


[0.635, 1.382] 


[0.616, 1.335] 


s 2 


[0.955, 1.072] 


[0.959, 1.075] 


[0.959, 1.072] 


[0.955, 1.075] 


[0.956, 1.072] 


[0.957, 1.072] 


% 

(derived) 


[0.376, 2.332] 


[0.546, 2.649] 


[0.554, 2.526] 


[0.416, 2.238] 


[0.452, 2.229] 


[0.503, 2.329] 


F0M3 


18.95 


22.94 


25.61 


5.989 


1.098 


6.755 


^3 


[0.439, 2.073] 


[0.509, 1.588] 


[0.544, 1.663] 


[0.093, 2.771] 


[0.122, 2.654] 


[0.159, 2.519] 


£3 


[0.851, 1.321] 


[0.708, 1.268] 


[0.818, 1.319] 


[0.903, 1.824] 


[0.546, 5.391] 


[0.871, 1.825] 


^3 

(derived) 


[0.069, 4.035] 


[0.129, 2.889] 


[0.244, 3.066] 


[0.052, 26.31] 


[0.051, 35.69] 


[0.037, 14.78] 


F0M4 


43.94 


62.03 


63.93 


8.878 


2.014 


9.214 




[0.553, 1.754] 


[0.588, 1.404] 


[0.593, 1.427] 


[0.061, 1.269] 


[0.032, 1.288] 


[0.052, 1.194] 


s 4 


[0.887, 1.148] 


[0.858, 1.134] 


[0.878, 1.144] 


[0.173, 1.552] 


[0.381, 5.181] 


[0.167, 1.563] 


(derived) 


[0.152, 2.708] 


[0.376, 2.424] 


[0.398, 2.458] 


[0.052, 18.46] 


[0.024, 169.9] 


[0.031, 22.27] 



TABLE IV: 95% C.L. for the parameters fj,i, rji, and Si from the second parametrization 57, 68] for redshift and scale dependence 
as a binned 2x2 parametrization with < z < 1, 1 < z < 2 and scale separation at k x — 0.01 hMpc' 1 and k x = 0.1 hMpc~ x . 
The data used is described in the text above and indicated as: ISW-galaxy cross correlations (ISW), weak-lensing tomography 
(WL). In all cases, the data used is combined with WMAP7 temperature and polarization spectrum (CMB), the matter power 
spectrum (MPK), SN, BAO, BBN, AGE and H . The FoM for each data set and {fn, £;} bin is also given. 

time). There is a slight decrease of the FoM for the {/Xi, Si} bin because there is not enough WL data in this bin to 
increase the FoM when added to ISW. Some tension between ISW and WL data sets is more explicitly seen in bins 
{(11, Si} and {(i2, S2} (see Table IV|) for both scale separators. Future surveys and other data sets may work better 
with scale separator at k x =0.1 hMpc~ 1 (45| . 

We give the results for the third binned method in Table IVl and show the relevant ellipses for the FoM's in Fig. [6] 
The constraints in Table IVl are tighter than the other binning results in Table ITVl and have higher FoM's corresponding 
to smaller ellipses in Fig. [5] The improvements in the constraints and FoM's comes from larger redshift bins allowing 
more data to constrain the parameters (called an accumulation effect by 57]) and these larger redshift bins also 
allowing a longer constant value for late times. The {(12, S2} bin is slighter looser in this third binned method 
because more of the ISW data points are in the {(11, Si} bin, now, due to a larger redshift bin, i.e. < z < 1.5, 
than in the previous two binned methods. We still observe some small tension when adding WL to the ISW data set 
combination with the FoM's and contour plots, and it is more easily seen for example in the {(ii, Si} and {fJ-2, S2} 
bins in Table IVl where the FoM's rather decrease. 

Before summarizing our results further below it is worth discussing the possible physical origins of the tensions 
on best-fit modified gravity parameters prefered by different combinations of data sets. To further explore this we 
decided to fix the core cosmological parameters to their WMAP7 values and use each set individually and removing 
the CMB data set since its data points will most likely dominate parameter constraints, and likely wash out some 
of the tension in our best-fit parameter spaces. Our results expectedly gave us weak constraints and for both the 
first parametrization and the second parametrization (in its binned form), we continued to observe tensions where 
they had previously existed (the tension in the unbinned parametrization was removed by fixing the core cosmological 
parameters). The best-fit parameter spaces obtained did give us further insight into the preferred parameter values of 
each data set. For the first parametrization we saw that the best-fit for MPK preferred the lowest value of s out of all 
the best-fits. Additionaly, in contrast to WL best-fits, it had negative best-fit values for R and V, though the best-fit 
Q and D parameters are actually not that different between. Combining the two data sets however we find that a 
best-fit for s that is higher than either of the two data sets individually, additionally a higher Q and correspondingly 
more negative V and R are preferred. For the binned second parametrization we looked for tension between the 
MPK+ISW and WL data sets. In the first bin, the WL data set had a much higher best-fit (i and S values than the 
MPK+ISW data set. The MPK+ISW+WL data set though has the lowest best-fit Si and a similar besti-fit (i x to 
the MPK+ISW data set again showing the tension when combining these data sets. These markedly different best-fit 
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values in the modified gravity parameter space between the individual and combined data sets highlight the tension 
but do not explain it. 

Working again with the first parametrization to explore the tension there, we next looked at the ouptut matter 
power spectrum for each of the best-fit parameter sets. This showed just how much an exacerbation the parameter s 
may be causing. The amplitude of the matter power spectrum produced with the best-fit WL parameter set was 2 
orders of magnitude greater than the amplitude of the spectrum produced by the best-fit MPK parameter set. Upon 
seeing this marked difference in amplitudes one can recall that in calculating the likelihood for MPK the amplitude of 
the power-spectrum is treated as a nuisance parameter [60| . This may be one reason for the tension when combining 
data sets, as a number of MG parameter combinations can produce matter-power spectra with similar shapes in the 
k-range where SDSS-DR7 data is present but with different amplitudes. The matter power spectrum is of course 
influenced by the growth of structure within a given gravity model. 

We recall again from above that the ISW-cffect is dependent on the quantity (4>+tp) 2 , while the weak-lensing signal is 
dependent on the quantity <p+ip. Furthermore, there has been much discussion in previous works [H], [H, Ht| as to how 
different values of the various MG parameters in the various parametrizations affect cosmological observations. We 
will recall some of that discussion related to the growth of structure in order to examine how different MG parameters 
will affect the matter-power spectrum. In the matter dominated era, the growth of structure (matter overdensities) 
on subhorizon scales (but still within the linear regime) is governed by the usual equation A m + HA m + k 2 ip = 0. 
Superhorizon growth of matter overdensities is a bit more complicated. We refer the reader to Eq. (A5) of (55j 
for a full expression of the growth equation in the matter dominated era (fj, there should not be confused with the 
ji in this paper as they are different). As mentioned in the previous reference a k 2 ip term will be present at all 
scales (it is actually scale- independent through Eq. (fTTT) '). This term is affected directly by the parameters fi or V 
from the second and third parametrizations respectively, and correspondingly the product of Q and R from the first 
parametrization. An increase (decrease) in these quantites from their GR values will increase (decrease) the magnitude 
of the gravitational potential tp. This increase (decrease) in the magnitude of the potential in turn combats (enhances) 
the suppression of the growth caused by the accelerated expansion of the universe, therebye enhancing (suppressing) 
the growth relative to its GR value on all scales. This will in turn increase (decrease) the overall amplitude of the 
matter power spectrum and the weak-lensing correlations relative to their GR values. Though not explicitly shown 
in the aforementioned equation for superhorizon growth, the next dominant term in the growth equation on those 
scales involves the derivatives of the potentials. Thus they will involve the derivatives of the MG parameters (these 
are shown in the mentioned equation though). In the second parametrization the derivatives do not matter much as 
they are zero for most of the evolution. For the first parametrization, however, the derivatives directly involve the 
parameter s as a multiplicative factor so the parameter s affects the growth on these scales in more ways than just 
govering how long the gravity modification is effective. 

For the second parametrization the effects of the parameters on the matter power spectrum is fairly straight forward. 
As discussed above an increase (decrease) in [i will increase( decrease) the overall amplitude of the matter power spec- 
trum. The parameter E however only affects the low-fc matter power spectrum. A S > 1(E < 1) will decrease(increase) 
the low-fc spectra relative to its E = 1 value for a given fi. This reflects the low-fc supprcssion(cnhanccmcnt) of the 
growth for those parameter values. In the binned form of this parametrization the parameters within a given fc-bin 
have the same effects as described above on the spectrum in that fc-bin, and do not affect the other fc-bin except near 
the transition between the bins. This allows parameters in different z-bins but the same fc-bin to be degenerate with 
one another in their effect on the matter power spectrum today. On the other hand, the effect of the first parametriza- 
tion on the matter power spectrum is a bit more complicated due to the presence of the parameter s. As discussed 
above, just as [i increases or decreases the amplitude of the matter power spectrum in the second parametrization, 
the value of the parameter combination QR will increase or decrease the amplitude of the matter power spectrum. 
How much the amplitude is changed though is affected by the value of s. A larger s means that the modification to 
gravity has not affected growth for very long, and thus the amplitude will be affected minimally. A smaller s on the 
other hand means the modification has been present for a much longer time and thus the amplitude of the matter 
power spectrum will be affected most for a smaller s. A smaller s will also affect the overall shape of the spectrum, 
mainly because when s < 1 gets smaller it allows the modification to affect the evolution of the potential outside 
of the matter dominated era. For the low-fc spectrum as with the second parametrization, a D > 1 (the parameter 
corresponding to E) will tend to cause a suppression in the low-fc spectrum relative to its D = 1 value (the converse 
still holds for D < 1. A larger s, however, tends to prefer a suppression in the low-fc spectrum and can even overpower 
the boosting effect of D < 1. 

Independent of parametrization changes to the matter power spectrum will of course affect the other observations 
such as the weak-lensing correlations, and the ISW-galaxy cross correlations. As discussed in [54j . for example, 
values of MG parameters moderately higher than their GR values can lead to an anticorrelation in ISW-galaxy cross 
correlations at low redshift. Those values of MG parameters promote growth while suppressing the ISW effect. A 
more detailed review of the effect of different MG parameters on the ISW-effect in the CMB TT power spectrum is 
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FIG. 5: 68% and 95% C.L. for the parameters and Si from the second parametrization [571 l6a ] for redshift and scale 
dependence as a binned 2x2 parametrization with < 2 < 1, 1 < 2 < 2 and scale separation at k x = O.lhMpc -1 . All 
contours are fit for CMB, MPK, SN, BAO, BBN, AGE and Ho- Contours enclosed by a solid line include WL, contours 
enclosed by a dashed line include ISW, and contours enclosed by a dash-dotted line include WL and ISW. The point (1, 1) 
indicates the GR values. 



discussed in [54j, [55[ . Basically both smaller and sufficiently larger values of MG parameters compared to their GR 
values will boost the ISW effect, while values of MG parameters moderately larger than their GR values will tend to 
suppress it. This is because the ISW effect is dependent on the square of <j> + ip. A smaller (larger) value for the MG 
parameters will increase (decrease) 4> + ip with larger values even causing the quantity to go negative. For sufficiently 
large MG parameters though <j) + tp will become so negative that the ISW effect is again boosted. 

The tensions for both parametrizations are possibly due to the way the MG parameters affect the shape and 
amplitude of the matter power spectrum and the way this change affects the other observables. 

In summary, our results using figures of merit and constraints from current data sets indicate that: 

(i) Using the first parametrization with functional form, the strongest current constraints on the MG parameters 
are found to be from the CMB+ISW+WL combination followed by the CMB+WL or CMB+MPK+ISW. 

(ii) We observe, in the case of the first parametrization that the combinations CMB+WL or CMB+MPK do better 
than CMB+MPK+WL. Similarly, the combinations CMB+ISW+WL or CMB+MPK+ISW do better than the 
combination CMB+MPK+ISW+WL. A closer look at the best-fit models preferred by each combination of data 
sets shows that tensions appear between the best-fit values for the MG parameters preferred by different data 
sets and that these tensions are much more pronounced when using this parametrization. 

(hi) Using the second parametrization or the binning methods, we find that the combination CMB+MPK+ISW of 
current data sets consistently provides some of the strongest constraints on MG parameters. 

(iv) When using the second parametrization or the binning methods, we find that the combination 
CMB+MPK+ISW+WL provides, at best, only little improvement when compared to the combination 
CMB+MPK+ISW, again indicating some tension between the data sets. But we find here that the tension 
between the best-fit MG parameters is much less pronounced than in the first parametrization although non- 
negligible. 
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Constraints for {fii, Si} binned parametrization < z < 1.5, 1.5 < z < 3 




^ = 0.01 


data sets 


WL 


ISW 


ISW, WL 


FoMi 


94.98 


120.9 


116.9 


Mi 


[0.629, 1.433] 


[0.537, 1.197] 


[0.540, 1.222] 


Si 


[0.929, 1.127] 


[0.948, 1.137] 


[0.939, 1.128] 


Vi 

(derived) 


[0.4004, 2.348] 


[0.702, 2.992] 


[0.659, 2.914] 


FoM 2 


149.6 


165.7 


163.5 




[0.620, 1.455] 


[0.570, 1.352] 


[0.579, 1.384] 


s 2 


[0.952, 1.071] 


[0.959, 1.074] 


[0.956, 1.069] 


V2 

(derived) 


[0.368, 2.307] 


[0.485, 2.616] 


[0.433, 2.542] 


F0M3 


28.20 


40.89 


46.09 




[0.393, 1.807] 


[0.594, 1.446] 


[0.604, 1.443] 


S3 


[0.869, 1.234] 


[0.795, 1.209] 


[0.862, 1.238] 


(d6riv6d) 


rn 1 q*3 a a 7oi 
[U.loo, 


rn /i Qoi 

JU.Zoo, Z.4oyj 


[U.oo/, Z.O/DJ 


F0M4 


48.17 


68.50 


69.97 




[0.494, 1.600] 


[0.545, 1.270] 


[0.574, 1.318] 


s 4 


[0.886, 1.145] 


[0.858, 1.133] 


[0.879, 1.146] 


(derived) 


[0.239, 3.135] 


[0.509, 2.687] 


[0.483, 2.631] 



TABLE V: 95% C.L. for the parameters Hi, r)i, and Si from the second parametrization [5^ . l68j | for redshift and scale dependence 
as a binned 2x2 parametrization with < z < 1.5, 1.5 < z < 3 and scale separation at k x = 0.01 hMpc -1 . The data used 
is described in the text above and indicated as: ISW-galaxy cross correlations (ISW), weak-lensing tomography (WL). In all 
cases, the data used is combined with WMAP7 temperature and polarization spectrum (CMB), the matter power spectrum 
(MPK), SN, BAO, BBN, AGE and H . The FoM for each data set and {m, Si} bin is also given. 



(v) In order to investigate these tensions, we fixed the core cosmological parameters and left the MG parameters 
to vary. In the cases of the second parametrization and the binning methods, the FoM's increase consistently 
(but only moderately in some cases) when adding a data set, indicating that some of the tension was reduced. 
A closer look at the best-fit MG parameters preferred by different data sets shows that, while reduced, a non- 
negligible scatter is still present. This explains the only moderate improvements obtained in some cases. This 
also indicates that there are non-negligible tensions intrinsic to the MG parameters preferred by different data 
sets. Moreover, using the first parametrization and fixing the core cosmological parameters did not reduce the 
strong tensions observed between the MG parameters which seems to suggest that the tensions are exacerbated 
by the functional form that the first parametrization is imposing on the MG parameters. 

(vi) We find that combining current data sets does not improve consistently the uncertainties on the MG parameters 
due to tensions between the best-fit MG parameters preferred by different data sets. Some functional forms 
imposed by the parametrizations can lead to an exacerbation of these tensions and this point requires further 
future investigations. It remains of interest to use these data sets separately in order to derive independent 
constraints on the MG parameters and thus allow one to cross- validate the results. 

(vii) Unlike previous work that used a binned parametrization and the CFHTLS lensing data, we do not find a 
deviation from GR using a similar method but the refined HST-COSMOS data, confirming thus previous claims 
in those studies that the previous result were likely due to some systematic effect in using CFHTLS lensing data. 

(viii) Finally, for all the parametrizations and binning methods used here, we find that the parameter values of general 
relativity are within the 95% confidence level contours for all data set combinations. 
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dependence as a binned 2x2 parametrization with < z < 1.5, 1.5 < z < 3 and scale separation at fc x = 0.01 h Mpc -1 . 
All contours are fit for CMB, MPK, SN, BAO, BBN, AGE and Ho. Contours enclosed by a solid line include WL, contours 
enclosed by a dashed line include ISW, and contours enclosed by a dash-dotted line include WL and ISW. The point (1, 1) 
indicates the GR values. 
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